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Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were 
conducted on the surface of a uranium tailing pond to better understand the differences between the two methods. 
An abnormal radon exhalation behavior was observed, leading to computational fluid dynamics (CFD)-based 
simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered. 
Based on the in-situ experimental and numerical simulation results, variations in the radon exhalation rate sub- 
ject to permeability, flow rate, and insertion depth were quantified and analyzed. The in-situ radon exhalation 
rates measured using the flow-through method were higher than those measured using the closed-loop method, 
which could be explained by the negative pressure difference between the inside and outside of the chamber 
during the measurements. The consistency of the variations in the radon exhalation rate between the experi- 
ments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of 
transient radon exhalation rates for diffusion and convection at the porous medium-—air interface. The synergis- 
tic effects of the three factors (insertion depth, flow rate, and permeability) on the negative pressure difference 
and measured exhalation rate were quantified, and multivariate regression models were established, with posi- 
tive correlations in most cases; the exhalation rate decreased with increasing insertion depth at a permeability of 
1x 1071! m?. CFD-based simulations can provide theoretical guidance for improving the flow-through method 


and thus achieve accurate measurements. 
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I. INTRODUCTION 


In the development of uranium resources or associated 
mineral resources containing natural radioactive isotopes, 
large amounts of solid radioactive waste are piled in tailing 
ponds or used as building materials. Radon exhaled by these 
materials poses potential risks to humans; therefore, radon 
exhalation has been extensively studied [1—6]. The radon ex- 
halation rate has been measured under various scenarios, such 
as soils, underground spaces, uranium tailing ponds, building 
materials, and indoor environments, to assess the potential 
risk of radon to the environment and humans [7, 8]. Accurate 
and reliable methods for measuring the radon exhalation rate 
are important. Currently, measurement methods for the radon 
exhalation rate include active (for example, closed-loop [9- 
12] and flow-through methods [13, 14]) ) and passive (solid- 
state nuclear track detectors [15, 16] and activated charcoal 
adsorption methods [17, 18]) methods. 

In the closed-loop method, the measurements obtained us- 
ing an accumulation chamber are influenced by the inser- 
tion depth, back-diffusion, and flow rate [19-21], whereas 
the negative pressure in the chamber plays a dominant role 
in the flow-through method [20-22]. The negative pressure 
difference between the outside and inside of the chamber is a 
byproduct of regulating the flow rate or deploying a radon pu- 
rification device at the chamber inlet. Generally, this factor is 


* This study was supported by the National Natural Science Foundation of 
China (Grant No. 11575080), the National Natural Science Foundation of 
Hunan Province, China (Grant No. 2022JJ30482), and the Hunan Provin- 
cial Innovation Foundation for Postgraduate (Grant No. QL20220206). 

t Corresponding author, yongjunye @ 163.com 


ignored in laboratory or in-situ experiments. For example, the 
higher radon exhalation rate measured via the flow-through 
method than via the closed-loop method can be attributed to 
measurement uncertainty, whereas the intrinsic difference be- 
tween the two methods has been ignored [14, 23]. Perme- 
ability characterizes the capability of a porous material to al- 
low fluids to pass through and varies greatly among differ- 
ent materials. Currently, the extent to which permeability af- 
fects measurements using the flow-through method remains 
unknown. Therefore, the effect of permeability should be 
quantitatively analyzed. 

Radon migration in porous media is affected by the radon 
diffusion coefficient, free radon production rate, permeabil- 
ity, and pressure gradient in a given medium. Although it 
is time-consuming and costly to regulate these parameters in 
laboratory or field experiments, the development of compu- 
tational fluid dynamics (CFD) techniques has allowed for the 
control of these parameters and the prediction of the trans- 
port of pollutants and flow field distributions in the soil, at- 
mosphere, and indoor environments [24—26], as validated by 
experimental data. 

In this study, small-scale measurements of the radon exha- 
lation rate on the surface of a uranium tailing pond were per- 
formed using the flow-through and closed-loop methods, and 
the differences between these two methods were assessed. To 
better understand the abnormally high radon exhalation rate 
measured with the flow-through method, CFD-based tech- 
niques were employed to simulate the radon-concentration 
field and dynamic transport in a porous medium and accumu- 
lation chamber, using the established three-dimensional ge- 
ometry and reasonable assumptions. Based on the numerical 
results, the effects of the flow rate, chamber insertion depth, 
and permeability of the porous medium on the transient radon 


exhalation rates for diffusion and convection were compre- 
hensively analyzed. The simulations may provide a theoret- 
ical basis for the optimization of flow-through measurement 
systems and devices to obtain accurate radon exhalation rates. 


Il. MATERIALS AND METHODS 
A. Radon migration theory 


The one-dimensional steady-state convective and diffusive 
migrations of radon in porous media can be described as fol- 
lows: 
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z Wr ei (1) 
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where C is the radon concentration, Bq m73, x is the ver- 


tical depth of the porous medium, m, D is the diffusion coef- 
ficient of radon in the porous medium, m? s71, v is the super- 
ficial velocity in the porous medium, m s7 t}, 77 is the porosity, 
%, a is the free radon production rate, Bq m~? s~!, and À is 
the decay constant of radon-222, that is, 2.1 x 107 s71 [27]. 
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Fig. 1. One-dimensional model of radon diffusion in porous media. 


The boundary conditions in Eq. (1) are C(0) = Ca and 
ƏC /ðx| —p, = Oat the top and bottom of the porous medium, 
respectively (Fig. 1). The radon concentration in the bound- 
ary air is significantly lower than that in the porous medium; 
therefore, Ca = 0 Bqm~® can be assumed. The analytical 
solution to Eq. (1) is 
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The radon exhalation rate can then be theoretically deduced 
as follows: 
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where Ea—e is the initial radon exhalation rate, which 
comprises diffusion and convection, Bq m?s—!, and C, is 
the radon concentration at the porous medium-—air interface, 
Bqm7?. 

For v = 0 ms !, Eq—e becomes Ea, which is the initial 
radon exhalation rate for diffusion because of the assumption 
of no convection inside the porous medium. 
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B. In-situ measurement parameters 


The radon exhalation rate was measured in the field on the 
surface of a uranium tailing pond located in southern China. 
A 1-kg sample of uranium mill tailings was collected within 
the 0-0.4 m depth range and mixed to obtain the following 
particle size distributions: <0.15 mm (10%), 0.15-0.3 mm 
(15%), 0.3-0.45 mm (24%), 0.45—-1 mm (49%), and >1 mm 
(2%). After sample drying, the measured particle density 
was 2500-+150 kg m~? based on the drained weight method. 
The mass radon production rate was 7.18 + 0.26 x 1074 
Bqkg~' s71, as determined by laboratory tests [27, 28] under 
dry conditions. 

A schematic of the flow-through measurement process is 
shown in Fig. 2. The inlet of the accumulation chamber was 
connected to the atmosphere via a vent tube. The outlet of the 
chamber was connected to a flow meter, desiccants, a radon 
detector (RAD7 from Durridge Company Inc., USA), and a 
pump used to regulate the flow rate. 

In the experiment, two cylindrical chambers with an inner 
diameter of @20 cm and a wall thickness of 5 mm were pre- 
pared. The heights of the two chambers were 11 and 13 cm, 
with insertion depths of 1 and 3 cm, respectively, suggesting 
that the effective height inside both chambers was 10 cm. The 
effect of the flow rate (0.5/1/2 L min~') on the radon exhala- 
tion rate was also considered. The RAD7 radon detector was 
calibrated at the Radon Laboratory of the University of South 
China prior to the experiment. 

A digital flow meter (MF4008 of Siargo Ltd., USA) and 
differential pressure meter (RS-YC-4G-4, Shandong Renke 
Control Technology Co., Ltd., China) were used to record 
the flow rate at the outlet of the chamber and the pressure 
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Fig. 2. (a) In-situ flow-through measurement device and (b)system 
connection diagram. 


difference between the inside and outside of the chamber, re- 
spectively. The soil temperature was 27+1.5 °C and the soil 
moisture was 12+2%. The soil moisture was monitored using 
soil sensors (Jinan Keyu Electronic Information Technology 
Co.). The RAD7 device was operated in the SNIFF mode 
with a cycle of 5 min for the closed-loop method and 3 min 
for the flow-through method. Sufficient time (approximately 
30 min) was reserved for purification of the RAD7 device af- 
ter each measurement. 


The accuracy of the digital flow meter was +1.5% of the 
full reading, with a resolution of 0.01 SLPM. The deviation 
between the two flow meters at the same flow rate was less 
than 1%. The accuracy of the digital differential pressure me- 
ter was +3% of the full reading +0.08 Pa at 25 °C, with a 
resolution of 0.01 Pa. The accuracy of the soil sensors was 
0.5 °C for the temperature reading and +2% for the mois- 
ture reading. 


The increase in the radon concentration in the accumula- 
tion chamber can be described as follows: 
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where t is the accumulation time, s, Æ is the radon exhala- 
tion rate, Bg m~?s~!, Sa is the inner surface area covered by 
the chamber, m?, Qout and Qin are the flow rates flowing out 
of and into the accumulation chamber, respectively, m? sl, 
V is the volume of the chamber, m3, , is the back-diffusion 
rate of radon in the chamber, s7}, and A; is the leakage rate 
of radon in the chamber, s~!. Moreover, Ae is the effective 
decay constant, which is equal to the sum of the radon decay, 
back diffusion, and leakage rates (\-=A+Ap,+)) s71. 

The boundary condition in Eq. (6) is C(0) = Co. The 
solution of Eq. (6) is 
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Once the radon concentration inside the accumulation 
chamber was stabilized, Eq. (7) became 
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By replacing A.V and QinCo with 0 as an approximation, 
based on the assumptions of E'S, > QinCo and Qout > 
AeV [14, 21], Eq. (8) becomes 
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Eq. (7) can be simplified to a linear form when (Qout/ V + 
Ae)t is sufficiently low. 
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For Co = 0 Bqm~°, the linear fit solution in Eq. (6) is as 
follows: 
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The linear duration of the increase in the radon concentra- 
tion for the flow-through method was 1 min owing to the high 
gas exchange rate in the chamber (a volume of 3.1416 L for 
the chamber, and a flow rate between 0.5 and 2 L min™ t). 

The exponential fit of Eq. (7) and the linear fit of Eq. (11) 
were applied during subsequent data processing to obtain the 
fitted radon exhalation rates of Ef—exp and Ef—lin, respec- 
tively, with the criterion R? > 0.9. R? depends on the resid- 
ual sum of squares (RSS) and total sum of squares (TSS), 
which can be calculated as 
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where n is the total sum of the data points, y; is the actual 
data point, f; is the value obtained from the fitted curve, and 
y is the mean value of all data points. 

The uncertainty of the radon exhalation rate can be deter- 
mined as [29] 
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where o E is the uncertainty of the radon exhalation rate, Æ 
is calculated as E = (A x B x V — QinCo)/Sa, and A (de- 
fined as os) and B (defined as Dout tV ), and their 
corresponding uncertainties 0 A and o B, respectively, can be 
obtained via exponential fitting. Because the measurement 
errors of V, Sa, and Qin were sufficiently low to be ignored, 
their contributions to o E were not considered. 


HI. NUMERICAL SIMULATION METHODS 
A. Geometric model 


In the numerical simulation, the cylindrical accumulation 
chamber had a wall thickness of 5 mm, and the internal effec- 
tive space of the chamber had a height and radius of 10 cm, as 
shown in Fig. 3. The two vent tubes (G6 mm) at the top of the 
chamber had a length of 8 cm, with four openings (Ø2 mm) 
in the walls of the tubes. Adjacent openings were spaced 2 
cm apart and distributed at 90°, and the resultant confined jet 
improved the uniformity of the radon concentration inside the 
chamber. The dimensions of the radium-containing medium 
were | m (length) x 1 m (width) x 0.5 m (height). 

To discretize the geometry, Fluent meshing was used to 
generate hybrid meshes that included polyhedral and hexa- 
hedral elements. The mesh was refined in regions where the 
radon concentration and gas flow changed rapidly (for exam- 
ple, the porous medium—air interface and vent tube), as shown 
in Fig. 4. 


B. Governing equations 


The gas flows in the chamber and porous medium were 
assumed to be incompressible and followed the laws of 
mass and momentum conservation. According to the law of 
mass conservation, the mass-conservation equation can be ex- 
pressed as follows: 


V-u=0, (14) 


where @ is the physical velocity vector, m s71. 


Following the law of momentum conservation, the 
momentum-conservation equation can be expressed as fol- 
lows [30]: 
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where pa is the density of air, kg m~°, P is the static pres- 
sure, Pa, g is the gravitational vector, m s-2, and F is the 
source term of the external body forces, Nm’, which de- 
scribes the viscosity of the porous medium and applies only 
to the porous zone. The laminar flow can be described as 
follows: 
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Fig. 3. Three-dimensional geometric model of the accumulation 
chamber: (a) Top view, (b) lateral view, and (c) overall view. 


where K is the permeability, m?, and p is the dynamic 
viscosity, Pa s. The viscous stress tensor (7) is described as 


[31] 


T= (u+ m) (V+ Va"), a7) 
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where u+ is the turbulent viscosity, m?s~!, which is de- 


Fig. 4. (a) Sectional plane of the discrete mesh where vent tubes are included. (b) Discrete mesh (overall view). 


fined as [31] 


k2 
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where C, denotes a constant, k denotes the turbulent ki- 
netic energy, m? s7?, and e is the turbulent dissipation rate, 
m?s~%, The governing equation for radon migration can be 
expressed as 
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Ot ae a dx; °° n’ ne 


nu 


where D is the radon diffusion coefficient, m? s7}, u; de- 
notes the physical velocity of convection along the x, y, and 
z directions, ms~', Spn denotes the source term, Bq m~? 
s71, and P; denotes the gas pressure along the x, y, and z 
directions, Pa. 

In a given porous medium, Spn is determined by the source 
and decay terms (Srn = —AnCrn + a), and D = Dm 
(Dm is the molecular diffusion coefficient of radon in air 
at 1.05 x 107° m? s71). In the accumulation chamber, the 


porosity was set to 1, Sr, was determined by the decay term 
(SRn = —AnCRn), and D = Dm + Dz (D: is the turbu- 


lent viscosity, m?s~!, which is defined as > H where the 


turbulent Schmidt number (Sc+) is 0.7 [32). 


C. Simulation parameters 


The simulated scenarios involved controlling the inser- 
tion depth (H1/H3/HS5 corresponding to insertion depths of 
1/3/5 cm), flow rate (0.5/1/2 Lmin~'), and permeability 
(1 x 1071/1 x 107-1°/1 x 107°? m?). The free radon pro- 
duction rate (œ) can be obtained as [27] 


a= ApsARaERn; (21) 


where a is the radon production rate, Bq m -¢ 4, Ps is 


the density of the porous medium, kg m73, Apra is the activ- 
ity concentration of radium, Bq ket, and Fp, is the radon 
emission coefficient, which is dimensionless. 

The radium activity concentrations in uranium mill tailings 
commonly range from less than 5 kBq kg to as high as 10 
kBqkg ` [33], with an emission coefficient ranging from 0.1 
to 0.35 [27], a dry bulk density of 1800 + 100 kg m~°, and 
an average porosity of 0.4, as measured for the samples in the 
laboratory. 


The diffusion coefficient is defined as [27] 


D=TDm, (22) 


where 7 is the tortuosity factor, 0.66 (dimensionless) [27, 
34]. 

The calculated diffusion coefficient of radon was 6.93 x 
107 m?s~!. Laminar flow in a porous medium follows 
Darcy’s law, and permeability, which is an intrinsic property 
of porous media, can be defined using the Kozeny—Carman 
empirical equation [35]: 
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where D, is the average grain size, mm, and K is the per- 
meability, m?. 

After sieving the tailing sand in the laboratory, D, was 
determined as 0.5 mm on average. K was calculated as 
2.47 x 107" me, 

Based on theoretical analysis and references, the deter- 
mined bulk density was 1800 + 100 kg m~3, the radium ac- 
tivity concentration was 7500 Bq ke, and the emission co- 
efficient was 0.1. Considering the linear relationship between 
the radon exhalation rate and free radon production rate in Eq. 
(5), the determined radon production rate was 3 Bg m~? s71 
in the simulations. Owing to the short duration (30 min) of 
the flow-through measurement and the low saturated radon 
concentration in the chamber contributing to a smaller back- 
diffusion effect than that of the closed-loop method, the de- 
termined radon diffusion coefficient was 6.93 x 1076 m? s71. 
Because permeability is a key parameter representing the 
ability of gas to flow through a porous medium, the perme- 
ability ranges from 1 x 1071! m? to 1 x 107°? m?, and the 
porosity is 0.4 based on references. The physical parameters 
for the simulations were ultimately determined, as listed in 
Table 1. 


Table 1. Physical parameters used in the simulation. 


Physical parameters Value 
Free radon production 3 (36, 37] 
rate, Bq m? g! 
Porosity 0.4 [28] 


Diffusion coefficient 6.93 x 107° [36] 
for radon, m? s7! 
Permeability, m? 1x10 yi x0 

/1 x 107° [38] 
Flow rate, L min™ t 0.5/1/2 


Insertion depth, cm 1/3/5 


The finite volume method (FVM) was adopted for the nu- 
merical calculations. The k — € turbulence model was used 
for the chamber zone owing to the confined jet from the four 
openings on the surface of the vent tube, whereas the laminar 
model was adopted for the porous medium. The empirical 
constants in the k — € turbulence model equations [31] were 


as follows: C, = 0.09, Cie = 1.44, Coe = 1.92, o% = 1.00, 
and o; = 1.30. 

To simplify the model, the following assumptions were 
made: the porous medium was considered homogeneous, and 
the impacts of gravity and temperature were ignored. 

For Qout © AeV and t > 4V/Qout (t>1500 s for Qout = 
0.5 L min™!) in Eq. (7), e~ '@eu/Y <1.9% indicates that the 
radon concentration in the chamber remained almost stable. 
A maximum duration of 1800 s was used in the subsequent 
simulations for the flow-through measurements. 


D. Initial and boundary conditions 


The initial radon concentration distribution in the porous 
medium was calculated prior to simulating the measurements, 
where the radon concentration at the porous medium-air in- 
terface was set to 0 Bqm~?. The measurement was con- 
ducted in transient mode with an initial radon concentration 
of 0 Bqm~® at 0 s in the chamber. The detailed boundary 


settings are presented in Table 2. 


Table 2. Boundary settings in the simulation. 


Boundary 
Inlet Pressure-Inlet 


Settings 
Radon concentration 
equals 0 Bq m?:0 Pa 
Equivalent to 0.5/1/2 
L/min, normal to 
boundary; radon 
concentration gradient 
equals 0 Bq m-4 
The surface of porous Pressure-Inlet Radon concentration 
medium not covered equals to 0 Bq m~; 0 
by the chamber Pa 


Outlet Velocity-Inlet 


Openings of the vent Interior — 
tube inside the 
chamber 
gas-solid interface Interior — 


Wall of the geometry Wall Radon flux equals 0 


Bqm 2s + 


E. Mesh independence 


Coarse, medium, and refined (0.5 M, 1.2 M, and 2.3 M 
cells, respectively) hybrid meshes were generated for the ge- 
ometric model. Subsequently, a mesh sensitivity analysis was 
conducted, as shown in Fig. 5(a). The radon exhalation rate 
approximately stabilized as the number of cells increased. 
The criterion for mesh discretization depends on minimizing 
the error between the numerical and theoretical results while 
reducing the calculation time. Thus, a medium mesh was used 
in subsequent simulations. 

The Reynolds number ranged from 89 to 517 according to 
the minimum (0.66 m s71) and maximum (3.8 ms~*) veloci- 
ties at the interface of the opening in the wall of the vent tube. 
By refining the mesh, the mesh sizes for the inlet, outlet, and 


openings of the tubes ranged from 0.12 to 0.3 mm, and yt 
reached a maximum value of 3.6 for the tube wall and 2.5 for 
the chamber boundary, which ensured that rapid flow near the 
wall was captured. 

The simulations were conducted using Ansys Fluent. The 
steady-state radon-concentration fields of the analytical and 
numerical solutions were compared before the transient sim- 
ulations. The analytical and numerical radon concentrations 
at different depths were determined, and the difference was 
found to be less than 0.2%. A detailed comparison is pre- 
sented in Fig. 5(b). 

The convergence condition was a residual of less than 
10~1° and a difference in the radon exhalation rate between 
two consecutive iterations of less than 1077. 
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Fig. 5. (a) Mesh sensitivity analysis. (b) Comparison of the radon concentrations at different depths between the analytical and numerical 


solutions. 


IV. DATA PROCESSING 


The volume-averaged radon concentration (C’,—ay,) in the 
chamber was obtained as follows: 


a (24) 


where n is the total number of cells inside the chamber, V; 
is the volume of cell i, mè, and C; is the radon concentration 
in cell i, Bgm7?. 

The volume-averaged static pressure (P,) inside the cham- 
ber was calculated as 


ae ee (25) 


where P; is the static pressure in cell V;, Pa. P, serves as 
the pressure difference under the 0-Pa condition at the porous 
medium—air boundary (outside the chamber). 

In the transient mode, the transient radon exhalation rate 
for the diffusive component can be defined as 


En—-a=nDVCle=0- (26) 


The transient radon exhalation rate for the convective com- 
ponent can be defined as 


En-ce=nuC]|x=0, (27) 


where u is the physical velocity for gas convection at the 
porous medium-air interface, m s™+. 

The overall transient radon exhalation rate is defined as fol- 
lows: 


En = En-a + En—e- (28) 


The radon exhalation rate on the surface of the porous 
medium before the measurement was obtained was referred 
to as the initial radon exhalation rate (Ea). For the homoge- 
neous porous medium in the numerical simulation, the initial 
radon exhalation rate was equal to the theoretical value ob- 
tained using Eq. (5). 


V. RESULTS AND DISCUSSION 
A. In-situ measurement of the radon exhalation rate 


Radon accumulation curves for the closed-loop and flow- 
through methods are shown in Fig. 6(a) and Fig. 6(b), re- 
spectively. The radon exhalation rate was exponentially fit- 
ted using Eq. (7), and R? varied between 0.92 and 0.96, as 


shown in Fig. 6(c). The radon concentration in the ambient 
air was monitored in all measurements, reaching 86 Bq m~? 
on average. The flow rates at the chamber inlet ranged from 
3.5 x 1076 to 1.7 x 1075 ms™t, as monitored using a digital 
flow meter. 

The radon exhalation rate measured using the flow-through 
method was higher than that measured using the closed-loop 
method and was positively correlated with the negative pres- 
sure difference and flow rate, that is, the higher the flow 
rate, the greater the radon exhalation rate or negative pres- 
sure difference, as shown in Fig. 6(d). An increase of 
161.88—293.2% (H3) was observed in the fitted radon exhala- 
tion rate determined using the flow-through method relative to 
that of the closed-loop method, with a slight increase at high 
flow rates between 1 and 2 Lmin~'. The negative pressure 
in the chamber was affected by the variation in the flow rate, 
with the negative pressure difference increasing by 433.98% 
for H1 and 671.43% for H3 from 0.5 to 2 L min™t. 


B. Numerical simulation of the measurements via the 
flow-through method 


1. Curve and linear fitting method 


To comprehensively examine the impacts of the flow rate, 
insertion depth, and permeability of the porous medium 
on radon exhalation during the flow-through measurements, 
CFD-based techniques were employed to simulate the flow- 
through measurements at insertion depths of 1/3/5 cm, per- 
meabilities of 1 x 10711/1 x 1019/1 x 107°? m?, and flow 
rates of 0.5/1/2 Lmin~!. The radon exhalation rate was ob- 
tained by applying exponential and linear fits, excluding the 
first 60 s of the data. The time ranges for exponential fitting 
were 60—600 s, 720 s, 900 s, 1020 s, 1200 s, 1500 s, and 1800 
s, whereas those for linear fitting were 60-240 s, 360 s, 480 
s, 600 s, 720 s, 840 s, 960 s, 1080 s, 1200 s, 1320 s, 1440 s, 
and 1560 s. A minimum duration of 60—600 s was considered 
for the exponential fitting because of the potential failure to 
fit data points with a shorter duration. The minimum duration 
of 60-240 s for linear fitting ensured that three-minute data 
points could be fitted after excluding the initial first minute of 
the data, as shown in Fig. 7. 

When a linear fit was applied, a short period of linear 
growth was observed in the flow-through measurements, as 
shown in Fig. 7. The slope of the linear fit at the same flow 
rate rapidly decreased with increasing time, whereas the slope 
differed for different flow rates. Therefore, a linear fit for the 
radon exhalation rate was not applicable to the flow-through 
measurements. 

In addition, the radon accumulation in the chamber at 1800 
s was close to saturation, and the steady-state radon con- 
centrations were 8457.508, 8500.607, and 9714.73 Bq ini 
at flow rates of 0.5, 1, and 2 Lmin“!, respectively. The 
radon exhalation rates were calculated to be 3.93, 7.06, and 
16.04 Bqgm~?s~! using Eq. (9) and 1.89, 2.94, and 5.79 
Bqm~?s~! using Eq. (8), respectively. The former val- 
ues were significantly higher than the latter, suggesting that 
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Fig. 6. Growth curves of the radon concentration for (a) the closed-loop method and (b) the flow-through method. Correlations between (c) 
the flow rate and fitted radon exhalation rate and between (d) the flow rate and negative pressure difference. 


AeV in Eq. (8) cannot be ignored in the flow-through mea- 
surements because this term, rather than the negative value in 
the closed-loop method, can reach a positive value in view of 
the observed negative pressure difference in the field experi- 
ments. Therefore, calculating the radon exhalation rate using 
Eq. (9) is also inapplicable to flow-through measurements. 

At t > 4V/(Qout + AeV) in Eq. (7), 1 — 
e~ (QourtrAcVItH/V > 0.981 for Qout = 0.5 Lmin™? indi- 
cates that the radon accumulation in the chamber had almost 
reached saturation. Therefore, a duration of 60-1800 s for 
the radon exhalation rate curve fitting was considered in the 
subsequent investigation. 


2. Effects of the permeability, insertion depth, and flow rate on 
radon exhalation 


To investigate the effects of permeability, insertion depth, 
and flow rate on radon exhalation in the flow-through mea- 


surements, the fitted (exponential and linear fits), transient, 
and initial radon exhalation rates were compared for different 
time ranges. Both the insertion depth and flow rate affected 
the radon exhalation; however, the flow rate had a greater ef- 
fect than the insertion depth, as shown in Fig. 8. With an 
increase in the flow rate from 0.5 to 2 L min™ t at an insertion 
depth of 1 cm and a permeability of 1 x 1071? m?, the curve- 
fitted radon exhalation rate (E f—exp) increased by 305.68%, 
whereas it increased by 427.66% for an insertion depth of 5 
cm. With an increase in the insertion depth from | to 5 cm 
at a flow rate of 0.5 L min™t, E’s—exp increased by 97.99%, 
whereas it increased by 137.09% at a flow rate of 2 Lmin™?. 
Both the increased flow rate and insertion depth contributed 
to the emission of radon at high concentrations, which con- 
tributed to an increase in the radon exhalation rate. 

Similar trends were observed for the variation in perme- 
ability, as shown in Fig. 9. With an increase in the flow 
rate from 0.5 to 2 Lmin~! at an insertion depth of 1 cm 
and a permeability of 1 x 1071! m?, Ef—ezp increased by 
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Fig. 7. Growth curves for the volume-averaged radon concentra- 
tion in the chamber (1 x 107~'° m?, H1). The slope of the linear fit 
varies with increasing flow rate, whereas the radon exhalation rate of 
porous materials should remain constant in a specific environment. 


135.66%, whereas it increased by 345.49% at a permeability 
of 1 x 107? m?. By increasing the permeability from 1 x 
1071 to 1 x 107° m? at a flow rate of 0.5 Lmin™, Ey_—cap 
increased by 251.70%, whereas it increased by 640.99% at a 
flow rate of 2 L min~'. The variation in the radon exhalation 
rate influenced by permeability could be conveniently simu- 
lated based on the assumptions of fixed porosity and radon 
diffusion coefficient. Changes in the permeability affect the 
ability of a porous medium to allow fluids to pass through. 
Higher permeability indicates lower fluid resistance, which 
promotes convective radon exhalation. 

Both En-a and En-—c gradually reached saturation after 
a certain period. In the case of a relatively low permeability 
(0-7 10-"m*); En-a was higher than En—c in most cases 
because the low permeability of the porous medium limited 
the radon at deeper layers from being drawn to the surface, 
and diffusive radon exhalation played a dominant role. In 
contrast, the case with higher permeability 1 x 107° m? was 
dominated by convective radon exhalation. Moreover, the in- 
creased flow rate had a significant effect on En-a and Eņn—c. 
Increasing the flow rate from 0.5 to 2 Lmin~* with a per- 
meability of 1 x 107!° m? resulted in increases of 353.26% 
for En—a and 871.34% for En—c at H1, and 434.65% and 
1243.71% for En-a and En—c, respectively, at H5. The larger 
the insertion depth of the chamber, the greater the impact of 
the increased flow rate on En—c than on En—a, where the high 
radon concentration around the insertion position is drawn to 
the surface and is responsible for the high convective radon 
exhalation rate. 

For an insertion depth of 1 cm, a permeability of 1 x 
10-7 /1 x 107° m?, and a flow rate of 0.5 Lmin™', the 
transient radon exhalation rate for diffusion (En—a) first de- 
creased and then gradually stabilized over time. This is be- 
cause radon accumulation in the chamber led to a decreased 
radon concentration gradient at the porous medium—air inter- 
face and inhibited diffusive radon exhalation. In contrast, the 
transient radon exhalation rate for convection (En—c) grad- 
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ually increased over time and reached equilibrium because 
of the higher radon concentration at the porous medium —air 
interface induced by convective radon exhalation. As En—c 
was not negligible relative to En-a, the transient radon exha- 
lation rate (En) was higher than the initial radon exhalation 
rate (Ea), and the difference increased with increasing per- 
meability, flow rate, and insertion depth. 

To quantify the effects of the permeability, insertion depth, 
and flow rate on the radon exhalation rate, En—e and Ep at 
1800 s, E f—exp for a duration of 1800 s, and F’¢_j;,, for a du- 
ration of 240 s were extracted, as shown in Fig. 10. The effect 
of varying the permeability from 1 x 1071! to 1 x 107° m? 
on En—e was significant, with increases of approximately 
3750% and 6100% at 2 and 0.5 Lmin™', respectively. At 
an insertion depth of 1 cm, the permeability increased by ap- 
proximately 251—650% for E f—exp and E'y_—jjn, as shown in 
Fig. 10a. A flow rate ranging from 0.5 to 2 Lmin™t in- 
creased En—c by approximately 870%, whereas it increased 
by 75-340% for Ef—exp and Ef—tin, as Shown in Fig. 10b. 
At 1 Lmin“', an insertion depth ranging from 1 to 5 cm re- 
sulted in a 58—162% increase in the radon exhalation rate, as 
shown in Fig. 10c. 

Ey,—c increased with flow rate or permeability at a fixed 
insertion depth, as shown in Figs. 8 and 9. Excessive con- 
vective radon exhalation during flow-through measurements 
was mainly responsible for the considerably higher measured 
radon exhalation rates than the initial values, which reason- 
ably explains the abnormally high radon exhalation rate in 
the field. 


C. Correlation between E, AP, and Q 


The previous line graphs revealed that both F,,_. and 
Ef—exp Were greater than Ea under variations in permeabil- 
ity, flow rate, and insertion depth. Considering the negative 
pressure difference in the chamber observed in the field ex- 
periments, streamlines were obtained from the simulations, 
as shown in Fig. 11 (an insertion depth of 1 cm, permeability 
of 1 x 10~!° m?, and flow rate of 1 L min™t at 1800s). The 
negative pressure difference in the chamber caused the gas 
outside the chamber to flow into the porous medium before 
entering the chamber, during which a high concentration of 
radon was transported into the chamber. This explains why 
E f—exp Was greater than Ea. 

Furthermore, the correlations between the flow rate and 
negative pressure in the chamber and E f—exp were obtained 
based on simulations, as shown in Figs. 12(a) and 12(b), re- 
spectively, indicating positive relationships. At a high flow 
rate of 2 Lmin!, the negative pressure difference (AP) for 
a permeability of 1 x 107°? m? increased by 779.4% (H5) and 
1202.67% (H1) when compared to that at a permeability of 
1 x 1071! m?. At a medium flow rate of 1 Lmin” t, AP in- 
creased by 500.16% (H5) and 767.31% (H1) but increased by 
332.92% (H5) and 495.74% (H1) at a low flow rate of 0.5 
Lmin~!. A comparable trend in the radon exhalation rate is 
shown in Fig. 12(b). At 2 L min !, E f—exp at a permeability 
of 1 x 107° m? increased by 640.99% (H1) and 1353.92% 


Le 
n 


a 
[e] 


— 
n 


Radon flux 
(Bq m?s’) 
& 


2 
n 


2 
= 


Time(s) 69 500 1000 1500 2000 
2.5 


id 
> 


— -00-90-9—9—9 


Radon flux 
(Bq m’s") 
> 


0.5 


0.0 
Time(s) 60 500 1000 1500 2000 


2:5 


ee 
S 


-— — — 
_ — 


Ur 


Radon flux 
(Bq m?s’) 
> 


S 
nv 


2 
S 


Time(s) 60 500 1000 1500 2000 
0.5 L min" 


60 500 1000 1500 2000 


60 500 1000 1500 2000 


60 500 1000 1500 2000 


12 


60 500 1000 1500 2000 
14.0 


60 500 1000 1500 2000 
14.0 


60 500 1000 1500 2000 
1L min! 2 L min! 


Fig. 8. Radon exhalation rates for the fitted values (E f—exp by curve fitting and Ef—ıin by linear fitting) under various durations, transient 
values in the numerical simulation (En component, where En-a and E,—- denote diffusion and convection, respectively), and analytical 
solution (Ea) (flow rate of 0.5/1/2 L min™ +, insertion depth of 1/3/5 cm, and permeability of 1 x 1071° m?). 


(H5) compared to that at a permeability of 1 x 1071! m?. It 
increased by 229.44% (H1) and 699.46% (H5) at 1 L mint, 
whereas it increased by 251.70% (H1) and 335.21% (H5) at 
0.5 Lmin~'. Similarly, when the flow rate increased from 
0.5 to 2 Lmin™t (H1), E f—ezp increased by 135.66% and 
345.49% at permeabilities of 1 x 1071! and 1 x 107° m?, re- 
spectively. Clearly, permeability, flow rate, and negative dif- 
ferential pressure profoundly influence the radon exhalation 
rate. 


The linear and nonlinear relationships between the nega- 
tive pressure difference and the variables, that is, the inser- 
tion depth, permeability, and flow rate, can be represented us- 
ing multivariate functions, as described by Eqs. (29) (R? = 
0.696) and (30) (R? = 0.978): 


Pmf = — 0.0541 x H 
+ 3.3943 x 109 x K, (29) 
— 3.1994 x Q 


Pf = —3.5591 x 1074 x 9.0787 x 7¢-0.3551 x 1-7366, 

(30) 

where P,,,¢ is the fitted value from the multivariate func- 

tions, Pa, H is the insertion depth, 1/3/5 cm, K is the perme- 

ability, 1 x 1071/1 x 1071/1 x 107°? m?, and Q is the flow 
rate, 0.5/1/2 Lmin™!. 

Similarly, the nonlinear relationships of E f—ezp with the 
insertion depth, permeability, and flow rate can be expressed 
using multivariate functions, as described by Eqs. (31) (R? = 
0.93) and (32) (R? = 0.983): 
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Fig. 9. Radon exhalation rates for the fitted values (E f—exp by curve fitting and Ef—ıin by linear fitting) under various durations, transient 
values in the numerical simulation (En component, where En-a and En-—e denote diffusion and convection, respectively), and analytical 
solution (Ea) (flow rate of 0.5/1/2 L min™ +, insertion depth of 1 cm, and permeability of 1 x 10 1 x 107+°/1 x 107°? m?). 


Ems =1.4632 + (1.959 x 10° x K x Q — 0.0761 x H) x H 
| 7.3055 x 108 x K 


+ 1.6846 x Q 
(31) 
Em = 10676 x Fy (1-21x 10° x K x Q—0.01x H) 
x K9:35579 , (32) 
7 @?-04 


where E,,,¢ is the fitted value from the multivariate func- 
tions, Bgm7?s71. 

In the simulation with a permeability of 1 x 1071 m?, 
E's exp decreased with increasing insertion depth, whereas 
for 1 x 107°? m?, the correlation was positive. Interestingly, 
in the case of 1 x 1071? m?, E f—exp initially increased with 


3 


increasing insertion depth before decreasing. This nuanced 
behavior was controlled by the terms (1.959 x 10° x K x Q— 
0.0761x H) in Eq. (31) and (1.21x 108 x K xQ—0.01x H) in 
Eq. (32). A goodness of fit approaching 1 for the multivariate 
nonlinear fit suggests that the three variables, insertion depth, 
permeability, and flow rate, synergistically affect the negative 
pressure difference and radon exhalation rate. Additionally, 
the radon exhalation rate could be predicted and corrected us- 
ing multivariate regressions. 

Both the simulation and experimental results suggest con- 
sistent trends in the variations in the radon exhalation rates 
and negative pressure differences under certain insertion 
depths and flow rates. Deeper insertion of the chamber re- 
sulted in a significant increase in the negative pressure differ- 
ence at a high flow rate of 2 Lmin™t. This amplification 
of the negative pressure difference subsequently increased 
the radon exhalation rate. Therefore, the numerical method 
demonstrated reliability in revealing trends in the radon exha- 
lation rates measured in the field. 
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Fig. 10. Growth rates of the transient radon exhalation rate (En), transient exhalation rate for convection (En-—c), fitted radon exhalation 
rate by curve fit (Ef—exp) and linear fit (Ef—in) subject to variations in (a) the permeability (from 1 x 107+ to 1 x 107° m? with a 1-cm 
insertion depth), (b) flow rate (from 0.5 to 2 L min~! with a 1-cm insertion depth), and (c) insertion depth (from 1 to 5 cm with 1 L min7*). 


Fig. 11. Streamlines inside the accumulation chamber and porous medium at 1800 s (1 x 107'° m?, 1 Lmin~', H1). Owing to the negative 
pressure difference inside the chamber, the gas outside the chamber passes through the porous medium before entering the chamber. In this 
process, a high concentration of radon from deep is transported into the chamber, resulting in a higher cumulative radon concentration in the 


chamber. 


D. Analysis of the variation in — Ae 


In Eq. (6), the term —A¢ can be interpreted as the radon 
dissipation inside the accumulation chamber, including radon 
decay, back-diffusion, and leakage. By applying an exponen- 
tial fit to Eq. (7), the value of Laut + Ae can be obtained 
to calculate the value of Ae; the results are listed in Table 
3. In previous analyses, certain relationships were observed 
between the negative pressure difference (or fitted radon ex- 
halation rate) and permeability, insertion depth, and flow rate. 
The impact of these three factors on radon exhalation is also 
reflected in the variation in Ae. 


In most numerical simulations (except for H1, which in- 
volved a permeability of 1071! m? and a flow rate of 0.5 
1 x 1071! m?), a positive value of —\, indicated the reverse 
process of radon dissipation, in which more radon flowed 
in and accumulated inside the chamber. In this process, the 
back-diffusion effect increases with the accumulation of high 
radon concentration in the chamber, reducing the increase in 
— àe. Specifically, by increasing the permeability from 1071 
to 107°? m? with an insertion depth of H5, Ae increased by 
413.13% at 0.5 L min™* and by 37.01% at 2 L min™+. When 
the flow rate increased from 0.5 to 2 L min~! with an inser- 
tion depth of H5, Ae increased by 2347.43% for 1071! m? 
and by 210.26% for 107°? m?. These data indicate that the 
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Fig. 12. Correlations between (a) the flow rate and negative pressure difference and (b) the flow rate and curve-fitted radon exhalation rate 
(Ef—ezxp, 60-1800 s). In the steady state before the measurement (that is, Q=0 L min‘), the pressure difference between the inside and 


2-1 


outside of the chamber is 0 Pa, and the initial radon exhalation rate is 1.46 Bgm ^ s 


term —A, had a positive value, and a higher permeability or 
flow rate caused a higher radon concentration in the chamber, 
as well as a nonnegligible back-diffusion effect at the porous 
medium-air interface, inhibiting the increase in Ae. 


Table 3. Term Ae calculated for certain cases of the insertion depth 
(H1/H3/HS), permeability (10~''/10~'°/10~°m?), and flow rate 
(0.5/1/2 Lmin“*). 

eS) 


: ILmin! 2Lmin 


Insertion depth Permeability(m?) 0.5 L min` 


H1 10 x 1071 3.97 x 10-°  —7.48 x 10`% —4.09 x 10-2 
10x 107'° 3.84 x 1074 —1.79 x 107? —4.66 x 107% 
10x 107°? -6.48 x 1074 —1.59 x 107 —3.40 x 107% 
H3 10x10-" —1.14 x 1074 —9.47 x 1074 —4.48 x 1073 
10x 107'° —5.42 x 1074 —1.97 x 107? —4.39 x 107° 
10x 107°? —7.61 x 1074 —1.48 x 1073 —2.53 x 1073 
H5 10x 1071! —2.14 x 1074 —1.13 x 107? —5.00 x 107% 
10 x 1071? —6.78 x 1074 —2.21 x 1073 —4.42 x 1073 
10x 107°? —8.84 x 1074 —1.46 x 107? —1.86 x 107° 


VI. CONCLUSION 


To explore the phenomenon of the considerably higher 
radon exhalation rate measured in field experiments using the 
flow-through method than the closed-loop method, a three- 
dimensional geometric model was established, and CFD tech- 
niques were employed to numerically simulate the radon 
transport and exhalation patterns in the accumulation cham- 
ber and porous medium under certain conditions of perme- 
ability, flow rate, and chamber insertion depth. The impacts 
of these factors on the radon exhalation rate were quantita- 
tively analyzed and compared. The conclusions are as fol- 
lows. 

The radon exhalation rate measured in the same chamber in 


the field using the flow-through method was higher than that 
measured using the closed-loop method. The difference in the 
values between the two methods increased with the flow rate, 
which can be explained by the increase in the negative pres- 
sure difference between the inside and outside of the chamber 
caused by the variation in the flow rate. 


The trend of the radon exhalation rate in the simulations, 
subject to variations in the flow rate and chamber insertion 
depth, was consistent with that of the field measurements. 
Using the simulations of the flow-through measurements, 
the evolution of the diffusive, convective, and total transient 
radon exhalation rates was obtained. The variation in the to- 
tal radon exhalation rate with the flow rate was closely related 
to the convective radon exhalation rate induced by the nega- 
tive pressure difference. In addition, applying a linear fit to 
the radon exhalation rates in flow-through measurements was 
not feasible because of the short time required for the radon 
concentration to stabilize. Additionally, the term —Ae in Eq. 
(6) could not be ignored when calculating the radon exhala- 
tion rate as a function of the saturated radon concentration 
during flow-through measurements, as Ae could reach a pos- 
itive value, which suggests that more radon flowed into the 
chamber because of the negative pressure difference in the 
direction opposite to that of leakage in the chamber during 
closed-loop measurements. 


Multivariate linear and nonlinear regression models were 
established for the fitted radon exhalation rate and negative 
pressure difference between the inside and outside of the 
chamber. Specifically, the negative pressure difference and 
fitted radon exhalation rates were positively correlated with 
the variation in permeability and flow rate, whereas the cor- 
relation with the insertion depth of the chamber depended on 
permeability. A negative correlation between the radon ex- 
halation rate and the insertion depth was observed at a per- 


meability of 1 x 1071! m?, whereas a positive correlation 
was observed for 1 x 107°? m?. At a medium permeability 
of 1 x 1071? m2, the radon exhalation rate first increased and 
then decreased with increasing insertion depth. The better 
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Nomenclature and units 


ARa Radium activity concentration, Bq kg~ 
G Radon concentration in the air, Bq m° 
Co Initial radon concentration, Bq m” 
Ca Radon concentration in the atmosphere, Bq m° 
Ch Radon concentration at the porous medium-air interface, Bq m`’ 
Cy—avg The volume-averaged radon concentration in the chamber, Bq m~”? 
D Radon diffusion coefficient in porous medium, ms! 
Dm Molecular diffusion coefficient of radon in air, 1.05 x 10-° m? s~! 
Dp Average grain size, mm 
Dt Effective turbulent diffusion coefficient, m? s~ 1 
E Radon exhalation rate, Bq m`? s7! 
Ea Initial radon exhalation rate equivalent to the analytical solution in diffusion, Bam ~? s~! 
Eag Initial radon exhalation rate equivalent to the analytical solution in diffusion and convection, Bq m 4s! 
Eş Fitted radon exhalation rate, Bq m'a! 
Ef—ezp Radon exhalation rate by exponential fit, Bq m -s ! 
Eflin Radon exhalation rate by linear fit, Bam ~? s~! 
Emf Radon exhalation rate obtained by the fitted multivariate function, Bq m?s! 
En Transient radon exhalation rate (En-a +En—-), Bq m 2s! 
En-a Transient radon exhalation rate for diffusion, Bq m -s ! 
Ente Transient radon exhalation rate for convection, Bq m 7s! 
Ern Radon emanation coefficient of porous medium, dimensionless 
F Source term of external body forces, N m` 
FVM Finite volume method 
H1, H3, and H5 | The insertion depth of the accumulation chamber into porous medium 
k Turbulent kinetic energy, m?s 7 
K Permeability, m? 
AP The pressure difference, Pa 
P Static pressure, Pa 
E; The volume-averaged static pressure in the accumulation chamber, Pa 
Pmt The fitted negative differential pressure, Pa 
Qout and Qin |Flow rate during measurements at the outlet and inlet, respectively, of the chamber, L min! 
Sa The area of the surface where radon exhales, m? 
Sct Schmidt number 
SRn Source term of governing equation of radon migration in numerical simulation, Bm ~? s~" 
t Time, s 
v Superficial velocity of gas flowing in porous medium, ms * 
V Volume of the accumulation chamber, më 
u Physical velocity of gas flowing in porous medium, ms * 
£ Vertical depth of porous medium, m 
Q Free radon production rate, Bq ms! 
E Turbulent dissipation rate, m?s’ 
n Porosity, % 
A Radon decay constant, 2.1 x 107" s71 
Àb Back-diffusion coefficient, s 
Ae Equal to A+Àb+A;, R 
AÀ The leakage rate for radon in the accumulation chamber, st 
H Dynamic viscosity, Pa s 
Let Turbulent viscosity, Pa s 
Ps Porous medium density, kg m~? 
Pa Air density, kg m” 
T Tortuosity factor, dimensionless 
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